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ABSTRACT 

As deeper observations discover increasingly distant galaxies, characterizing the properties of high-redshift 
galaxy populations will become increasingly challenging and paramount. We present a method for measuring 
the clustering bias of high-redshift galaxies from the field-to-field scatter in their number densities induced 
by cosmic variance. Multiple widely-separated fields are observed to provide a large number of statistically- 
independent samples of the high-redshift galaxy population. The expected Poisson uncertainty is removed from 
the measured dispersion in the distribution of galaxy number counts observed across these many fields, leaving, 
on average, only the contribution to the scatter expected from cosmic variance. With knowledge of the A-Cold 
Dark Matter power spectrum, the galaxy bias is then calculated from the measured cosmic variance. The results 
of cosmological N-body simulations can then be used to estimate the halo mass associated with the measured 
bias. We use Monte Carlo simulations to demonstrate that Hubble Space Telescope pure parallel programs 
will be able to determine galaxy bias at z > 6 using this method, complementing future measurements from 
correlation functions. 
Subject headings: surveys — methods: statistical — galaxies: statistics 



1. INTRODUCTION 

Galaxy bias, the spatial clustering strength of galaxies rel- 
ative to the matter field, encodes important information about 
the connection between the luminous baryons of galaxies and 
the dark matter (DM) halos that host them. The galaxy bias b, 
defined as the ratio between the galaxy and matter correlation 
functions (CFs) b^ = £,gg/£,mm, reflects both the clustering of 
galaxy-scale DM halos and how luminous galaxies populate 
halos of a given mass. A self-consistent connection between 
the luminosity function (LF), angular C F, and mass-to-light 
ratio of galaxies can be developed ("e.g.. Ivan den Bosc h et alJ 
2003: Vale & Ostriked 12004 IConroy et alJ i2006: Lee et alJ 
200l . 

CF estimates of bias, which require large samples 
of galaxies spanning a wide range of angular sepa- 
rations, are r emark a bly successful at l ow-redshift (e.g., 
iNorberg et all |2001|; IZehavi et alj |2002|) . Angular CF 
measurements on smaller samples have been performed 
for Lyman Break Galaxies (LB Gs^ (iGiavahsco et all 119981: 
Arnouts et al. 1999; Ouchi et all |2004 120051: iFoucaud et al " 
2(303; Hamana et al. 2004; Adelberger et al. 2005; Alle ri et al 
2005; Kashikawa et al. 2006; Lee et al. 2 006, 20095~and 
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Lyman-g emitters (lOuchi et al.l 120 03: Ha mana et alJ 

iKovac et al]|2007l; iGawiser et alj|2007i) over a range of red- 
shifts. These studies have extended high-significance (> 3a) 
detections of bias out to z ^ 5. 

At the highest redshifts (z > 5.5), limited samples 
have made s i gnifica nt detections of galaxy bias difficult. 
iMcLure et alJ (l2009l) used a 0.63deg^ area observed by the 
UKIRT Infrared Deep Sky Survey Ultra Deep Survey and the 
Subaru XMM-Newton Survey to measure the clustering of 
5 < z < 6 LBGs with z < 26A B (b = 5.2'^i), but mo st of their 
systems lie at redshifts z < 5.5. lOverzier et al.l (120061) reported 
a measurement of bias for /-dropout galaxies in Great Obser- 
vatories Origins Deep Survey (GOODS) (b = 4.i;2;g) but no 
detection of clustering in the Ultra Deep Field. A study of 
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z-d ropout galax ies in the Subaru Deep Field and GOODS-N 
by lOuchi etal] (f2009) indicated that z ~ 7 galaxies may be 
strongly clustered, but measuring bias from only 22 candi- 
dates is challenging. 

Given the difficulty of measuring clustering of high-redshift 
galaxies from the angular CF, an exploration of comple- 
men tary approaches is wa rranted. As a salient exam- 
ple, Adelberger et al.l (Il998l) used a counts-in-cells analysis 
(Peebles 1980) to measure the clustering of redshift z ^ 3 
LBGs from statistical fluctuations in their number densities. 
Having obtained spectroscopic redshifts for A^ ^ 270 LBGs in 
six fields, they split their samples into redshift slices and cal- 
culated the variance in galaxy counts in rectangular volumes. 
The bias was then measured from the excess dispersion in the 
number counts beyond the Poisson variance. 

In this Letter, we present calculations that demonstrate the 
bias of high-redshift (z > 6) galaxies can also be measured 
from field-to-field dispersion in galaxy counts from photomet- 
ric dropout samples. Specifically, we use Monte Carlo simula- 
tions to demonstrate that pure parallel programs with the Hub- 
ble Space Telescope (HST) that obtain large numbers (> 100) 
of statistically-independent samples of z ~ 6-7 galaxies with 
the Advanced Camera for Surveys (ACS) or Wide Field Cam- 
era 3 (WFC3) (e.g., Ratnatunga 2002; Sparks 2002; Tren^ 
2008; Yan 2008) can provide a significant measurement of 
high-redshift galaxy bias. Throughout, we adopt the cosmo- 
logical parameter values determined by the 5 - year Wilkinson 
Microwave Anisotropy Probe joint analysis dKomatsu et alj 
2009). 

2. MEASURING HALO BIAS FROM COSMIC VARIANCE 

Consider a number A^fieids of telescope pointings on the sky. 
At the location x, of the /-th field the expected number density 
of galaxies is 

n(x/) = n[l+MDM(x/)], (1) 

where h is the mean average galaxy number density, b is the 
bias, and (5dm = iP~P)/p) is the matter overdensity of the 
field. The matter variance (S^f^), averaged over many survey 
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volumes V, is 



>jIm(V,z)={SI^)=d\z) 



(2^ 



p{k)\w{Ky)? (2) 



where D{z) is the growth fact or at redshift z, P(k) is th e matter 
power spectrum (we use the lEisenstein & Hul Il998l transfer 
function), and 

lV(k, y) = sine f^^j sine (^^j sine (-^j (3) 

is the Fourier transform of the rectangular volume V = rQ^ x 
rQy X 6r with angular area Q^ x 9, and depth Sr at distance 
r. The expected dispersion in galaxy counts owing to cos- 
mic variance is s imply g^ = b^a^M i^- R; Stark etal. 2007; 
iTrenti & Stiavellill2008l; OvTunoz et al.ll2010;.Robertsoa.20ia) . 
With a sample of A^fleids fields, the distribution of galaxy 
counts A^, with mean N per field, will have a dispersion a^ = 
{(N-N)^) w a^N^+N that includes contributions from both 
cosmic and Poisson variance. The bias can be estimated as 



-A^ 



" "dm 



(4) 



Below, we demonstrate that this estimate can be used to mea- 
sure the bias of high-redshift galaxies from multiple widely- 
separated HST pointings. The accuracy of this method will 
depend on A^seids and the relative size of the cosmic and 
Poisson variances, and we now explore these effects through 
Monte Carlo simulations. 

3. MONTE CARLO SIMULATIONS OF GALAXY BIAS 
MEASUREMENTS FROM COSMIC VARIANCE 

To simulate high-redshift galaxy number counts, we require 
a model for the distribution of galaxy counts for each field in 
our sample. Below, we describe our model distribution and 
our Monte Carlo realizations of this distribution that act as 
our model observational samples. 

The initial cosmological linear overdensity field is a Gaus- 
sian random field with dispersion a^f^ given by Equation |2l 
As the density field evolves owing to gravitation, the over- 
density distribution changes to maintain the positivity condi- 
tion on the density ((5dm > -!)■ The shape of the quasi-linear 
overdensity distribution can be approxim ated by a lognormal 
dColes&Jonesll 19911: iKofman et al][T994 ). The correspond- 
ing galaxy overdensity distribution will be lognormal, written 
as 
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where y = (1 + Sg) and x = (1 + <JZ) ' ■ The expected cosmic 
variance is Cg < 1 (see below), so the quasi-linear distribution 
is appropriate. 

At each location x,, the number of galaxies observed will be 
a Poisson-sampled random variate of Equation |5] The prob- 
ability of observing A^ galaxies in fields of volume V with 
cosmic variance cTg, given the expected number N = hV , is 

_ \ r°° _ - 

p(N\N,al)=— / d6,pi^^(S,\al)[(l+6,)Nfe-^'^'^''"' (6) 



dAdelberger et alJll998l) . Our simulations of the bias mea- 
surement consist of drawing discrete random samples from 
this distribution, as described below. 



3.L Observational Model 

Measuring the bias from the dispersion in galaxy counts re- 
quires an observational program suited to probing many high- 
redshift galaxy samples. Our presented method is motivated 
by HST programs (Programs 9488, PI Ratnatunga, 9575 and 
9584, PI Sparks, 11700, PI Trenti, and 11702, PI Yan) that 
use pure parallel observations with ACS or WFC3 to probe 
z > 6 - 7 galaxy populations. While the total numbers of 
long-duration and multiple-passband observations in the HST 
archive are few (~ lOi), we expect that similar pure paral- 
lel observations will be obtained by HST in ongoing and fu- 
ture cycles. We will therefore model future HST pure parallel 
programs that obtain ACS /- and z-band data, which would 
allow for /-dropout selection at z ^ 6 to z "^ 26.5A B (sin- 
gle orbit exposure, extrapolated from iBeckwith et an i2006). 
Additionally, as an example comparison with WFC3 obser- 
vations, we will model the Yan WFC3 program that can se- 
lect F600LP-F098M dropouts at 6.5 < z < 7.3, and reach 
F098M ^ 27 AB magnitude sensitivity (single orbit exposure), 
but our method can be adapted for use with any WFC3 filter 
strategy. We adopt these survey designs, summarized in Ta- 
ble[T] as our points of comparison for modeling observational 
measurements of high-redshift galaxy bias. 

We wish to simulate the number of high-redshift galaxies 
these pure parallel^ surveys will discover in each field. The 
mean abundance N per field is determined by integrating a 
ISchechten (Il976h LF over the volume probed by each point- 
ing. For modeling the ACS observations of /-dropouts at 
5.5 < z < 6.5, we adopt the measured z ^ 6 LF parameters 
0.0014 Mpc"^m ag"', M^ = -20.24AB, and a = -1.74 
:007h . With these values, we expect A^ = 4.3 
per 11.3 arcmin"^ ACS field in the magnitude range Muv < 
-20.5. To approximate the conve rsion between za b appar- 
ent and restframe Muv used by Bouwens et alJ (I2006I) . we use 
Muv-Zab~-47. 

For modeling the WFC3 observations the z '^ 7 LF is un- 
certain at 6.5 < z < 7.3, but if we adopt the measured z ^ 7 
LF parameters ^^^ = 0.0014 Mpc" ^mag"', M^ = -19.91AB, 
and a = -1.77 (lOesch et all 1201 Ol) . we expect A^ = 1.04 per 
4.7 arcmin^ WFC3 field in the UV absolute magnitude range 
Muv ^ -20.2. Our r esults wil l be simila r for other estirn ates 
of the z - 7 LF (e.g., !McLu re^et al. 2010i: lYan et ai1l2009h . 

Galaxy counts for each field can be generated by discretely 
sampling the Equation|6l with a variance a„ and bias b deter- 
m ined using t he abundance-matching methodology presented 
in iRobertsonI ( l2010i) . This method assumes the galaxies oc- 
cupy DM halos wit h an abundance and bia s appropriate for 
their high redshifts (iTinker et al.ll2008ll2010h . 

3.2. Uncorrelated ACS Fields 

The left panel of Figure [T] shows A^reais = 100,000 realiza- 
tions of the galaxy bias measurement for z ^ 6 galaxies with 
M, 



(IBouwens et 



UV 



<- 



-20.5AB across Nmds = 30-200 (light to dark blue) 
widely-separated ACS fields. For this example magnitude 
range at redshifts 5.5 < z < 6.5, the mean number of galax- 
ies per 11.3 arcmin^ is iV = 4.3 and ctdm ~ 0.05. The model 
bias calculated for this galaxy abundance is b ^ 1 .5 (verti- 
cal dashed line). The mean galaxy count per ACS field N 
and the sample variance a^j are calculated from each realiza- 
tion. The "observed" bias is then estimated using EquationHl 
histogrammed, and normalized such that the area under each 
curve is unity. 
For each A^geids-sized sample the mean expected measured 
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TABLE 1 
Model Parallel Surveys 



Redshift 



Detector Orbits/Field" 



Depth 
[AB mag] 



[AB mag] 



0-. l/VN Fields/SWR ~ 3*= 



5.5<z<6.5 
6.5<z<7.3 



ACS 
WFC3 



z ~ 26.5 
F~27 



1.4 X 10-3 d 
1.4 X 10-3 = 



-20.24'i -1.74'i 
-19.91"= -1.77"= 



4.29 
1.04 



0.37 
0.50 



0.48 
0.98 



40 
130 



" Includes single orbits for the detection, dropout, and bluer veto bands. 

Average number of galaxies per field. 
*= Number of fields required to reach SNR ~ 3. 
■^ Adopted from jBouwens et al.l 120071) . 
■= Adopted from lOesch et alj <201Ct) . 
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Fig. 1. — Monte Carlo simulations of galaxy bias measurements from the field-to-field vaiiation in number counts across N = 30-200 (light to dark blue) 
statistically-independent, widely-separated fields. Shown are 100,000 realizations of the measurement for z ~ 6 ACS /-band (left panel) and z ~ 7 WFC3 
F600LP-band (right panel) dropouts with UV absolute magnitudes Muv ^ — 20.5AB and Muv ^ — 20.2AB, respectively. As the number of independent fields 
increases, the SNR of the cosmic variance measurement increases rapidly. For so me re alizations the scatter in the distiibution of number counts is smaller than 
Poisson and cosmic variance cannot be measured (b = Q bins, for a discussion see j|3.2) . 

bias is accurate, with A/, = {{b)-b)/b sa -0.025 for A^fieiAs = 30 
decreasing to A^, w -0.014 for A^fieids = 200. The expected un- 
certainty ah in the measured bias decreases from cri,/b w 0.36 
for A^fieids = 30 to (Jb/b « 0. 16 for A^geids = 200, improving as 
oc l/VA^eids- When A^fleids is small, the measured dispersion 



a^^ can be less than A^, providing b^ <Qm Equation |4] Such 
catastrophic failures happen infrequently for ACS /-dropout 
samples (< 7.5% for A^fieids = 30 and < 0.3% for A^feids = 100, 
shown as the b = Q bin). 

Given these results, we suggest that using the dispersion in 
galaxy counts determined from ACS photometric samples in 
uncorrelated fields to measure bias will be a novel method for 
learning about the spatial clustering of z ^ 6 galaxies. The 
bias measurement significance increases from SNR ~ 3 with 
A^fieids = 40 to SNR w 6.4 with A^fcids = 200 fields for /-dropouts 
with magnitudes near M*, providing a method for measuring 
bias that increases in sensitivity with an increasing number of 
independent samples. 

3.3. Uncorrelated WFC3 Fields 

The right panel of Figure [T] shows A^reais = 100,000 realiza- 
tions of the bias measurement forz ^ 7 galaxies with absolute 
magnitudes M\ 



uv 



<- 



-20.2AB across A^fieids = 30-200 (fight to 
dark blue) widely-separated WFC3 fields. For this magnitude 
range at redshifts 6.5 < z < 7.3, the expected mean number 



of galaxies per 4.7 arcmin^ is A^ = 1 .04 and ctdm ~ 0.06. The 
model bias calculated for this galaxy abundance is /? w 8.5 
(vertical dashed line). The observed mean galaxy count per 
WFC3 field A^ and the sample variance cr^ are calculated from 
each of the A^reais realizations. The "observed" bias is then esti- 
mated using Equation m histogrammed, and normalized such 
that the area under each curve is unity. 

These histograms demonstrate that characterizing the bias 
of rare z ^ 7 galaxies will be challenging. Approximately 
Mieids =130 pointings are needed to measure the bias of these 
galaxies with an appreciable signal-to-noise {SNR ^ 3) in an 
accurate way ( A^ < -0.023). With fewer fields (A^fieids < 100), 
the expected uncertainty in the measured bias is large ((Jb/b > 
0.36) and overestimated (A^ > 0.05). When A^fieids is small, 
the measured dispersion a^ can be less than N, providing 
b^ <0 in Equation!!] In these catastrophic failures (26% for 
A^fieids = 30 and 9% for A^fieids = 100, shown as the /? = bin), 
an unusually small Poisson scatter limits sensitivity to the cos- 
mic variance in number counts. These failures decrease as the 
number of pointings increases (to 3% for A^geids = 200). 

The larger numbers of z '-^ 6 galaxies per ACS field com- 
pared with the number of z ^ 7 galaxies per WFC3 field (see 
^3.2b makes the z ^ 6 ACS experiment much easier. While 
Mieids ^130 pointings are required to significantly (SNR ^ 3) 
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measure the bias, such a determination would provide a cru- 
cial constraint on the nature of z ^ 7 galaxies. Current HST 
programs may have too few fields to perform this experiment, 
but extensions or new pure parallel programs will likely pro- 
vide enough pointings in the next few HST cycles to success- 
fully measure the bias of z ^ 7 from cosmic variance. 

3.4. Correlated Fields 

Pure parallel observations with ACS and WFC3 provide a 
novel method to obtain statistically-independent samples of 
the high-redshift galaxy population. However, the forthcom- 
ing HST Multi-Cycle Treasury (MCT) programs will probe 
greater numbers of contiguous pointings. While the MCT ar- 
eas may be used to measure the angular correlation function 
of high-redshift galaxies, one may wonder whether contigu- 
ous areas may be used to measure bias from cosmic variance. 
Here, we estimate the effects of correlations in the galaxy 
counts induced by large-scale density modes shared between 
nearby pointings. We consider z ^ 6 dropouts in ACS point- 
ings, as the improved statistics relative to z ^ 7 dropouts more 
clearly illustrate the effects of field-to-field correlations. 

The detailed correlation between nearby fields will depend 
on the geometry of the tiling and the separation between 
fields, and are difficult to calculate for general survey designs. 
However, we estimate the field-to-field correlatio ns in a con- 
tiguo us tiling as follows (for another estimate, see iHu & CohnI 
l20()6h . 

Consider as before A^fieids pointings, but arranged in a con- 
tiguous tiling. The cosmic variance of an independent field 
with volume V will depend on al(V,z)- Correspondingly, the 
galaxy counts in the larger volume A^fieidsV will have a sam- 
ple variance crg(A%eidsV,z). The number counts within the i-th 
field of the A^fieids pointings within the contiguous survey will 
then be correlated with the other j fields approximately at the 
level of 

p- — ~^~27Tr~^- ^^^ 

This typically underestimates the correlation between very 
nearby or adjacent fields, but should approximate the aver- 
age correlation of the most widely-spaced pointings within 
the contiguous area. The matrix describing the correlation be- 
tween fields / and j then has a simple form, with the diagonal 
elements E,, = 1 and off-diagonal elements S/jy, = p- 

To model the counts in correlated fields, we use a stan- 
dard Cholesky decomposition technique to instill the distri- 
bution of galaxy overdensities with the correlation matrix E, 
and then Poisson-sample the correlated distribution to gener- 
ate the counts in each field. Figure|2]shows the results of this 
procedure. Since the number counts in fields are correlated 
(p w 0.3 for A^fieids = 30 and p w 0.09 for N^d. = 200), the dis- 
persion a^ in the ACS galaxy counts becomes smaller than in 
the uncorrected case. The resulting estimates for the cosmic 
variance are correspondingly smaller, leading to a bias mea- 
surement that is underestimated (A;, = -0.15 for A^fieids = 30 
and Afc = -0.06 for A^feids = 200) and is more susceptible to 
catastrophic (measured b <0) failures (13.6% for A^eids = 30). 

3.5. Estimating Halo Masses from Galaxy Bias 

The simulations presented in ^3.2l and[ 3.3l demonstrate that 
the bias of high-redshift galaxies can be measured directly 
from the variance in number counts between uncorrected 
fields. Using the connection between bias a nd DM halo 
mass calculated from cosmological simulations (iTinker et al.l 
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Fig. 2. — Same as left panel of Figure 1, but for a contiguous survey where 
large scale structure correlates the number counts in each field. The field- 
to-field coiTelations decrease the dispersion in galaxy counts, and leads to a 
10-20% underestimate of the galaxy bias. 

I2OO8II2OIOI) . the galaxy bias can be used to estimate the char- 
acteristic mass of DM halos with similar spatial clustering 
strength. For detailed discussions of the c onnection between 
hal o mass and b i as, we refer the reader to iRobertsonl (1201 Oh 
and lTinkeretaP (l2010h . 

Figure [3] shows the characteristic hal o mass estimated by 
converting the bias (Figure [U using the [Tinker et al.l (1201 Oh 
bias-mass relation determined from cosmological simula- 
tions. For the case of Muv S -20. SAB magnitude /-dropout 
galaxies at redshift z ^ 6 with bias b « 7.5, the mass of ha- 
los with the same bias at z ^ 6 is Mhaio « 2.1 x 10"/i"'Mq 
(dashed line). The bias can provide a precise estimate of a 
characteristic DM halo mass (8% in logM for A^geids = 30 and 
3% in logM for A^geids = 200) that, while skewed, is accurate 
(Am = [(logM) -logM] /logM = -0.03 for A^fieWs = 30 and 
Am = -0.007 for A^geids = 200). This measure of a characteris- 
tic halo mass only corresponds to the halo mass of /-dropout 
galaxies if each halo hosts o ne galaxy and m ass strongly cor- 
relates with luminosity (see IRobertsonl 12010 ). but provides 
a convenient conceptual tool for discussing the approximate 
mass scale of high-redshift galaxies. 

4. SUMMARY AND DISCUSSION 

We have presented a simple coun ts-in-cells method (e.g., 
lPeebleslfT980l: lAdelberger et al.ll 1 9981) for measuring the bias 
of high-redshift galaxies from the field-to-field variation in 
number counts induced by cosmic variance. The number of 
high-redshift galaxies of a given luminosity are measured in 
a large number independent, widely-separated fields. The 
Poisson contribution to the variance in the number counts 
across these fields is removed, leaving a cosmic variance 
contribution that depends on the bias and the mean matter 
overdensity fluctuations on the scale of the survey field vol- 
ume. We use Monte Carlo simulations of the bias measure- 
ment to estimate the effectiveness of the method for example 
galaxy populations, /-dropouts at z ~ 6 in ACS-sized fields 
and F600LP-dropouts at z ^ 7 in WFC3-sized fields. Our 
method can be adapted for other dropout filter strategies. At 
z ^ 6 (z ^ 7) the method should provide a measurement of 
galaxy bias with SNR > 3 for luminosities near Muv ~ M^ 
when A^fieids > 40 (A^fieids ^ 130) fields are used. The uncer- 
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Fig. 3. — Characteristic halo mass infeired from the bias of z ~ 6 ACS 
r-band dropouts with UV absolute magnitudes M\]y < — 2().5AB. The bias is 
measured from number count variations across N = 3Q — 200 (light to dark 
blue) fields (see left panel of Fig[TJ. and then converted into a mass using the 
rKnker et al. 1 2010) bias function for DM halos. 

tainty on the galaxy bias improves with an increasing number 
of independent fields as « l/^/Nmds■ These requirements on 
the number of fields are concomitant with the expected num- 
ber of pointings availabl e from extensions to o ngoing pure 
pai-allel HST programs ( Trentill2008l; 1^120081) . and the po- 
tential for using this method to measure high-redshift galaxy 
bias is therefore promising. We also show that using corre- 
lated, nearby fields to perform the presented measurement 
typically leads to an underestimate of the bias owing to a 
correlation-induced decrease in the dispersion of galaxy num- 
ber counts. 

The measured bias can be associated with a characteristic 
mass of DM halos with similar clustering. We show that 



the typical uncertainty in the bias measured using the pre- 
sented method corresponds to an uncertainty in the inferred 
DM halo mass of ^ a few percent in logM for /-dropouts. 
Combining these estimates of halo mass with the observed 
luminosities allow for an estimate of a mass-to-light ratio, 
which can constrain the conne ction between the galaxies' ob- 
serve d LF and their clustering dYang et al.ll2003l : IConrov et alj 
|2006|) . Numerous effects can alter the correspondence be- 
tween bias and halo mass, and we therefore view the pre- 
sented method as primarily a measure of bias and interpret 
the halo mass estimates with caution. Our simulations assume 
that the high-redshift galaxy population is essentially volume 
limited and pure. Methods to handle incomple teness for LBG 
samp l es have already been e ngineered (e.g., iBouwens et alj 
120081: iReddv & Steide l"2009'). Contamination at a fractional 
level of fr will lower t he bias by an amount A^ « ( 1 - fc) 
(e.g.. lOuchi et al.|[2004l) . Using simulations that include red 
z ^ 2 galaxy interlo pers with fc = 0.2 and a bias of /? w 5.5 
(lOuadri et alJl200'7h . we find that the measured bias at z ^ 6 
is lowered by ^ 17%. This systematic error is comparable to 
the statistical error on the bias when A^fieids J^ 200. Lastly, not 
every parallel field will be useful for measuring high-redshift 
galaxy counts owing, e.g., to possible bright star or Galactic 
reddening contamination, and the efficiency of the presented 
method may be correspondingly decreased. 
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